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We study the B°-B° mixing amplitude in Standard Model by computing the relevant hadronic 
matrix element in the static limit of lattice HQET with the Neuberger light quark action. In the 
quenched approximation, and after matching to the MS scheme in QCD, we obtain B^(mi 7 ) = 
0.940(16) (22). 
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1. Introduction 

B® — B® mixing is highly important in testing the Standard Model (SM) and constrains strongly 
its extensions. Since it is a flavor changing neutral process it occurs through loops so that the 
corresponding mixing amplitude is a sensitive measure of \V ts \ and \V t b\, as the major SM loop 
contribution comes from ?-quark. The mixing of weak interaction eigenstates B® and B° induces a 
mass gap AM S between the mass eigenstates B s h and B s i. Experimentally, only a lower bound to 
AM S is currently known, namely AM S > 14.4 ps _1 [1], and the hope is that experimenters will soon 
provide us with an accurate measurement. 

Theoretically the — B® mixing is described by means of an Operator Product Expansion, i.e. the 
Standard Model Lagrangian J^sm is reduced to an effective Hamiltonian ^^ =2 , up to negligible 
terms of 0(1 /M%): 

^fi =2 = ^K(V,lV ts ) 2 VBSo(x t )C(Vb)Qf L =2 (Vb), Qf L =2 = b7nLsb^LS, Hb~m b (1.1) 

where Tjg = 0.55 ±0.01 So(x t ) is a known Inami-Lim function of x t = mj/M^, [2], C(Hb) is the 
Wilson coefficient computed perturbatively to NLO in a s (flb) in the MS (NDR) scheme, and <2l£ =2 
is a four-fermions operator coming from the reduction of the box diagrams in Jz?sm to a local oper- 
ator in the effective theory. The hadronic matrix element of g££ =1 is conventionally parameterized 

aS — R 

<*°ie£! =2 (M»K> = ^I/bBbM . (i-2) 

where B^Xl^b) is the B s meson bag parameter and fg s its decay constant. 

So far BeXfJ-b) nas been computed by using lattice QCD [3]-[9]. One of the major problems with 
those computations is in the following: the standard Wilson light quark lattice action breaks ex- 
plicitely the chiral symmetry, which tremendously complicates the renormalization procedure of 
Qll =2 an d its matching to the continuum. To get around that problem we compute BsXl^b) by 
using the lattice formulation of QCD in which the chiral symmetry is preserved at finite lattice 
spacing [10]. On the other hand, it should be stressed that our heavy quark is static, as the currently 
available lattices do not allow to work directly with the propagating b quark. Thus our results will 
suffer from 1/ra;, -corrections. 

2. Computation on the lattice 

In our numerical simulation we choose to work with the action S = Sjp 1 + Sf, where 

5 ™ = « 3 L{ /_J+ W [ h+ ( x ) -V HYP1 >-6)/i+(;t-6)] -h-(x) [V? YP (x)h-(x + 6) -h-(x)] } 

is the static limit of HQET action [11] which has been modified after using the so-called HYP 
(hypercubic blocking) procedure [12], that is enough to substantially improve the signal/noise ratio 
[13] [the field h + (h~ ) annihilates the static heavy quark (antiquark)]. Sf = a 3 J2 r \j/(x)D ( j!J l i//(x) is 
the overlap light quark action with 

£>JT ) = (l-^am )D N + m , D N = ^(l + ^L=j, X = D W -^, 
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Figure 1: Effective binding energy of the -state when currents are local (unfilled symbols) or smeared 
(filled symbols) 



where Dw is the standard Wilson-Dirac operator. The overlap Dirac operator D N verifies the 
Ginsparg-Wilson relation {y 5 ,D^ l °^} = ^D^^^D^ 10 ^ and the overlap action is invariant under 
the chiral light quark transformation [14] 

which is essential to prevent mixing of four-fermion operators of different chirality [15]. In other 
words, in the renormalization procedure, the subtraction of the spurious mixing with d = 6 opera- 
tors will not be needed. 

We thus compute the two- and three-point functions: 

££h0 = <E^o (^^(0)>„ ^2 A e~«, (2.1) 



C { y} +AA (ti,t) = (E^&fOOi&OjV&O)* ZAviBslO^nMh e- £ ^-'\ (2.2) 



f,— f»0 



,-e(ti-t) 



(2.3) 



x.y 



Av=h ± 'fiy i s, d 1 =h^%(l-y s )s i h^ j Y^-fy j , 2 = /* (+)i (l -yV#-)'(l -yV'. 

\/Z4 = (0|Aq \B s ) v = (0|Aq \B s ) v and e is the binding energy of the pseudoscalar heavy-light me- 
son. In C^ 2 '* (ti,t) one current Aq is local whereas the other is smeared. The role of the smearing 
is to isolate earlier the ground state [16], as shown in Fig. [I] 1 . We see that the same state is isolated 
when purely local currents are used (with those currents the signal does not exist if V HYP is not used 
in the heavy quark action). The source operators in C VV+AA (?;,?) and C ss+pp (tj,t) are the smeared 



currents Aq , whereas the four-fermion operators 0\ and 2 are purely local. In (|2~l|), (1.2) and (2.3) 



the subscript "v" and superscript "~" are designed to remind the reader that states and operators 
are denned in HQET. Note that in the computation of C vv+AA (ti,t) and C ss+pp (ti,t) there are two 
terms, coming from two different Wick contractions, namely ^Ba(t)^B jj(ti) and ^Bij{t)B ji(ti), 



where i, j are the color indices and = Tr £ f y^l^l 1 " (0;x,t)yoy 5 (x,f,0) 
are the light and heavy propagators respectively and the trace is over spinor indices. 



; y L and 2> E 



'Even if the time interval from which we extract the binding energy starts at t = 9 (green line), the overlap with 
radial excitations is quite reduced since t = 6 when currents are smeared. 
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Figure 2: Signals for R\<i{ti,t) defined in eq (2.4): green and blue lines indicate the time interval on which 
we fit the signal to a constant to extract By (a) and B^ (a) respectively 



After having computed the correlation functions (2.1), (2.2) and (2.3) we build the following two 
ratios Ri(tt,t) and R 2 (tut): 



R\(U,t) 



t ,-,»o v{B s \di\B s ) v 



l -Z 2 C 



Ri(ti,t) 



c 



(2)- 
A^AA 

(3) 



SS+PP 



(ti)C%t(t) 
{tut) 



fl(0|Ao|B,) v | 2 

ti-t>0 v(B s \Q2\B s ) v 
~* -f|(0|A |B. v ) v |2 



B { {a), 



= B 2 a). 



(2.4) 



Those ratios are calculated either with a fixed time t G [—6, — 8, — 10, — 12, — 14, — 16] and t\ free, 
or by fixing t\ G [6,8, 10, 12, 14, 16] while letting t free. We take the average of the two options. In 
Fig. H we show the quality of the signals for Ri^(tj,t), with t\ = 6 fixed. The signal for B~i(a) is 
quite stable as a function of t t , whereas the signal for B 2 {a) rapidly deteriorates for larger tj, and is 
completely lost for t\ > 10. 



3. Extraction of physical Bg s 

Three steps are required to extract Bb s =B\ from the lattice: 

(1) B\ 2{a) are matched onto the continuum MS(NDR) scheme at NLO in perturbation theory at 
the renormalization scale \i = \ja [15], 

(2) §12 are evolved from ;U = \ ja to jU = m\, by using the HQET anomalous dimension matrix, 
known to 2-loop accuracy in perturbation theory [7, 17], 

(3) #i,2 (M = m b) are then matched onto their QCD counterpart, B\^{mb), in the MS(NDR) scheme 
at NLO [17]. 

The advantage of using a chiral light quark action for the step (1) lies in the fact that four-fermion 
operators can mix only with a four-fermion operator of the same chirality. In other words we have 
not more than 4 independent renormalization constants in the renormalization matrix, because 0\ 
and 02 can mix neither with O3 = h^J^LS hry^s, nor with O4 = h + (l — y 5 )s h~( \ + y 5 )s : 

Bf\ll) \ = ( Z n (aH) Z 12 (an)\ ( B { {a)\ 
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p 


V 


N CO nf 


P 




6.0 


16 3 x 32 


80 


1.4 


0.06 GeV 


5.85 


16 3 x 32 


30 


1.6 


0.09 GeV 



Table 1: Parameters of our simulations: OTq and p have been chosen following [18, 19] 

Actually, because of the heavy quark symmetry, those constants are not all independent. By using 
the heavy quark symmetry (HQS) transformations h^(x) —^i^i^^j.^ (/ = 1,2,3), 

and the equations of motion for the heavy quark h^Jo = ±/jW, we see that: 

HQS(i) HQS(i) 
OSS+PP = —O[yV+AA) : Ow+AA ► Oyv+AA, Oss+PP ► — 0(yV+AA) t - 

As the action is invariant under the HQS transformations, we can deduce important constraints on 

the renormalization matrix Z/y: 

(0 V v+aa{h)) = Z n {O vv+AA (a)) +Z n (0 S s+pp(a)), 
(Ovv+aa(h)} = Z n (Ovv+AA(a)) -Z n (0 {vv+AA)i {a)) (HQS(i)), 

which implies that Z\i = 0. Moreover 

(Oss+pp{IJ-)) = Z 2 i(0 V v+AA(a)} +Z 2 2(0 S s+pp(a)), 
-<<Vv+AA),(A0) = Z 2l {O vv+AA {a)) -Z 22 {0 [vv+AA)i {a)) (HQS(i)), 

" £ {vv+aa)i (^)±0 {vv+aa)q (^i) = -(Oss+pp^)) - (O vv+AA (n)), 

i=l3 

= (3Z 2 i -Z 22 ){Ovv+AA{a)) -Z 22 (O ss+P p(a)}, 
= -(Zn +Z 2l )(O vv+A A(a)) -Z 22 (0 S s+pp(a)), 

which leads to Z 2 \ = (Z 22 — Z\\)/A. 
Finally we have 

Bp (At) \ = / Zn(afl) o 

^ S (A0 J \ 1 [Z22(an)-Z u {afl)] Z 22 {a\l) J \B 2 {a) J 

Therefore only two independent renormalization constants are required to match the bag parameters 
2*1,2 (a) computed on the lattice to their counterpart renormalized in MS scheme. 

4. Results and discussion 

Our results are based on two simulations, with the parameters given in Tab. |I[ We find B^ s (mb) = 
0.940(16) (22), where the first error is statistical, the second is systematic and contains the error 
from the estimation of a s (l/a) and the finite a effects. From Fig. || it can be seen that our value 
is larger than the previous static result [4]. This difference is likely due to the use of Neuberger 
light quark action (no subtractions), due to the use of the HYP procedure, or the combination of 
both. From Fig. || we also notice that our value is also somewhat larger than the results obtained 
with the propagating heavy quark, which is due to our neglect of 1 / my corrections or their not so 
proper renormalization. JLQCD collaboration showed that the errors due to quenching are likely to 
be small [8, 9]. We also plan to address that issue by unquenching the B° s — mixing amplitude in 
the static limit and by avoiding the subtraction procedure as well. The feasibility study by means 
of twisted mass QCD is underway. 
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Bernard, Blum and Soni (1998) 

Gimenez and Reyes (1999) 

APE (2000) 

UKQCD (2001) 

SPQcdR (2002) 

JLQCD, NRQCD (2003) 

JLQCD (unq. N f = 2), NRQCD (2003) 

Orsay, static heavy quark (2005) 



Figure 3: Various lattice values of Bg S {tnb) [3]-[9]; blue symbols correspond to a computation made with 
a static heavy quark 
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